Do Bose metals exist in Nature? 
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Abstract. We revisit the concept of superfluidity in bosonic lattice models in low dimensions. 
Then, by using numerical and analytical results obtained previously for equivalent spinless fermion 
models, we show that the gapless phase of ID interacting bosons may be either superfluid or - 
remarkably- metallic and not superfluid. The latter phase -the Bose metal- should be, according to 
the mentioned results, a robust and stable phase in ID. In higher dimensionalities we speculate on 
the possibility of a stable Bose metallic phase on the verge of a Mott transition. 
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INTRODUCTION 

In the last decades there have been a lot of numerical and theoretical works on interacting 
Bose gas in lattice or continuous models, fl 0, IH 0, H, El 0] The recent advance in 
the realization of optical lattices, where bosons are trapped on particular lattice sites 
and the interaction and the hopping parameters can be tuned continuously, has also 
opened a novel possibility to understand fundamental questions of many-body quantum 
mechanics, that can be experimentally checked with high degree of reliability and 
reproducibility. An important example is the realization of a Mott insulating state in 
a system with strong on site repulsion. [ 8, |9|] 

In this work we want to focus on an even more fundamental question, that is related to 
the concept of superfluidity. This concept deserves some discussion and generalization 
when considering a lattice model. The following discussion is not at all academic 
because, at present, lattice models can be realized with laser optical techniques, and 
the gedanken experiment we will discuss in the next section can be in principle realized 
experimentally. 

THE MODEL ON A RING 

We consider a one dimensional Bose-Hubbard model in the ring shown in Fig.dTJ). The 
lattice ring is rotating with given angular velocity COq with respect to the environment 
E which is considered here at rest for simplicity. Indeed in actual experiments the 
environment is usually rotating, but this does not change the forthcoming analysis, 
because our choice is just related to the reference frame. 
The Hamiltonian can be generally written as: 

Hv = Y, £ k4+ v a k+v+y (i) 
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FIGURE 1. The model: the site positions are rotating with respect to the environment with given 
velocity v. The rotating sites are decoupled from the environment E, whereas the bosons can interact 
weakly with it, but are allowed to occupy only the site positions of the rotating ring. 

where a^. = Y,Re^ k+v ^ R a R creates a boson in the ring with momentum (h = 1) k + v, 
£k is the dispersion of bosons in the lattice (e.g. £ k = —2tcosk for nearest neighbor 
hopping), periodic £k+2n = £k an d even e k = e_ k , V any two-body interaction term 
depending only on the relative distance between bosons [e.g. V = U /2Y,rHr(hr — 1) 
for the Boson-Hubbard model, where hr = ci\cir\, thus is unaffected by the velocity 
v = (OqL/k of the rotating frame. The total momentum in the reference system where 
the ring is at rest is given (modulo 2%) by: P = Y.k^ a l+ v a k+v an d the momenta k are 
obviously quantized according to the known relation kL = Inn. Strictly speaking in a 
lattice only the operator e ,p is defined, but this does not change the forthcoming analysis. 
In the forthcoming sections Hq will be indicated by H for simplicity. 

The experimental issue to detect superfluidity is related to the following experiment. 
After an experimentally accessible time (the ring rotating and the environment at rest) 
will all the bosons be at rest relative to the environment (or equivalently will they move 



with an appropriate velocity with respect to the ring lattice positions)? If this is not the 
case we can speak about supefluidity, a fraction p s of all bosons decouple from the rest 
and remains uncoupled from the environment. 

It is clear from the previous definition that superfluidity is related to the coupling 
to the environment (otherwise any finite momentum will be conserved for ever in the 
ring). Nevertheless it is possible to obtain a result that is independent of the interaction 
between the environment and the ring if the following three conditions are satisfied: 

i) the thermodynamic limit L — > °o is considered, 

ii) a finite temperature is given and the low temperature limit is considered after that 

the thermodynamic limit is employed, 

iii) the model Hamiltonian provides a stable phase in the low energy spectrum, namely 
stable for small perturbation of the Hamiltonian itself. 

The first two conditions are easily understood: only within the finite temperature 
canonical distribution Z = Tre~$ H the momentum e lP can equilibrate even without 
considering the coupling with the environment, and the probability of each eigenstate 
of the isolated ring H is given correctly by e~$ E[ , for a macroscopic system (L — > 
oo), just when the coupling environment-ring is negligible with respect to the bulk L. 
The coupling environment-ring is used only to equilibrate the system and obtain a 
property-superfluidity- that characterizes the system itself and not its coupling with the 
environment (otherwise we could talk about "superfluidity of capillary tubes" and not 
superfluidity of e.g. He<\). In order to achieve this consistent definition the Hamiltonian 
itself describing the system without environment has to define a stable phase of matter, 
namely a phase stable for small physical perturbations of the Hamiltonian, otherwise, 
clearly, the realization of a particular phase can obviously depend on the coupling 
environment-system. 

In cold atoms experiments L, the number of sites, can be as large as 10 5 and the 
thermodynamic limit at fixed temperature represents a realistic limit. It is important to 
emphasize that the physical zero temperature limit is highly non trivial in this respect. If 
we take first /3 — > °° and then L — > °° superfluidity cannot be tested because the lowest 
eigenstates of the Hamiltonian with non zero current have also a non trivial complex 
momentum e' p that is obviously conserved and no relaxation process can occur to the 
real ground state in a finite size system. If we explicitly consider a coupling system- 
environment as in fllCXI to induce current relaxation, it is clear that this process should 
be essentially equivalent to work in the thermodynamic limit with an arbitrary small 
temperature. 

We conclude therefore that the correct limit for detecting zero temperature superflu- 
idity is to take first L — > °° and then /3 — > °°. This order of the limits leads indeed to 
the definition of superfluidity that is independent of the coupling system-environment, 
whenever this is possible, namely when (iii) is satisfied. 



FREE ENERGY AND THERMODYNAMIC EQUILIBRIUM 



In this section we revisit some basic notions in thermodynamics, by introducing the basic 
quantities that define the superfluid density. The following considerations are completely 
general and hold in any dimensionality D with minor changes, that we omit in the 
following. 

In the thermodynamic equilibrium it is easy to show that the free energy: 

F v =-l/j81og(7>e-0 H ") (2) 

does not depend on v for discrete values of v = v n = Ikhl/L, with ni being an arbitrary 
integer. For large size L this limitation is very weak because we can reach any finite 
velocity v for L — > oo as , in this limit, the discrete velocity values merge in a contin- 
uum. Indeed for each L a simple unitary transformation, commuting with the two-body 
interaction V: 

U = e -iL R v n Rn R (3) 

removes the velocity v n from H. This follows immediately after simple application of 
canonical commutation rules, implying that [n R i,ciR\ = —8 R R /aR, so that U*ak +v U = a&, 

and finally U^H V U = H is easily obtained (£/"' 'a\ +v cik+ v U = (U^al +v U)(U^ak +v U) = 

a\cik in the kinetic energy). 

Using the above relation, it follows that the free energy: 

F v = -l/piog(Tre- pHv ) = -1/jB log(7>[/V^[/) = -l/piog(Tre- pu ^ HrU ) = F 

(4) 

does not depend on v whenever v = v n , namely when U is defined. Notice that in the first 
step we have used the invariance of the trace under cyclic permutation. 

From this relation we can expand the partition function in powers of v because the 
mapping k — > k + v n is just a shift of the finite size momenta and the kinetic energy 
of H v can be recasted in the following form: H = Y*k e k-v n a \ a k We thus obtain upon 
performing simple differentiations: 



dF 
~dv 



~<J> v + <K>v + O(v 2 )=0 (5) 

e(k 
dk 



= L^T a k a k ( 6 ) 

k 



^ d z e(k) t 

K = I^r^K ( 7 ) 



k 

where the brackets < O > v (< O >) denote the finite temperature averages 

_ TrOe^ Hv 
>v " Tre-f* H » 



on the Hamiltonian of the rotating ring (non rotating ring, i.e. with v = 0). Strictly 
speaking the previous differentiation in the free energy is not allowed because the 



possible velocities are quantized v n = 2nn/L otherwise the unitary transformation U 
is not properly defined. In a more rigorous way one can indeed consider that: 

d£ 

< J > v , =< tfjU > = < Y,-^\k+Vn4 a k > (8) 

k 

In the latter equation we can expand ^\k+v„ — ~w~ + v ^ dfi an ^ use mat < ^ >= ^ 
in the canonical ensemble of the rotating ring because the current is odd under reflection 
and the Hamiltonian H is even for v = 0. This immediately implies the linear relation 
between the current flowing in the frame rotating with the ring and the corresponding 
velocity at thermal equilibrium: 

<J> v =<K>v n + 0(v 2 n ) (9) 

within "weak" assumptions on the average boson occupation in momentum space = 

a\ak (e.g. 1 jL^k d ' jfy < n k > is finite) that allows to neglect the 0(v 2 ) term even for 
small but macroscopic velocities v n . It has to be remarked here that the fundamental 
relation © is valid only at thermal equilibrium and this may be obtained only after an 
exceedingly large time. This is indeed the case when, for L — ► oo, superfluidity occurs. 
On the other hand whenever the relation © is fulfilled the current flowing in the ring is 
just representing the condition of thermal equilibrium: all the bosons by scattering with 
the environment eventually converge to an equilibrium state characterized by no charge 
flow in the environment frame. 

We notice that a linear relation between the current and the velocity can be obtained 
within the linear response theory. The evaluation of < J > v for small v is given by: 



v (10) 

where J(t) = e tH Je~ tH . and can be obtained by simple expansion of the trace with simple 
and standard manipulations. Whenever the kernel Cf(0) relating the current response to 
an arbitrary small velocity is not equal to < K > we will have a relation current velocity 
plotted in Fig©. For any measurable finite velocity quantized as multiples of 2n/L, 
there is no net current flow in the environment frame, implying that < J > v — < K > 
v = 0(v 2 ) at equilibrium, as expected. However for unphysically small values of the 
current the linear response may have a finite slope as shown in Fig.©. 



< J > v = 



Jdt<J(t)J(0) > 



Dynamical limit co — » 

We are arguing in the following that the situation displayed in Fig.© is actually the 
common one for a superfluid (L^oo finite temperature). The point is that in a superfluid, 
in order to obtain the equilibrium steady state solution where no net current is flowing in 
the environment frame, an exceedingly large time is necessary because an initial current 
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FIGURE 2. Equilibrium current in the environment frame as a function of the velocity v of the rotating 
ring displayed in Fig.Q. The linear behavior occurs in an irrelevant region with exceedingly small 
velocity, certainly not measurable for large L, and in any case without any macroscopic current flowing in 
the environment frame. 



can very slowly relax to the steady state. In order to understand this fact it is useful 
to remark that the experiments of superfluidity are usually done with time dependent 
velocity v (e.g. a torsional pendulumflllll). If the frequency 0) is much larger than the 
inverse relaxation time of the current we can safely assume that linear response theory 
can be applied and we should obtain in this case a perfectly linear relation between the 
frequency dependent current and the frequency dependent velocity: 

J(cd) = a(co)v((o) (11) 

for either small or macroscopically measurable velocities as shown in the Fig.©. The 
situation is in some sense similar to the evaluation of the conductivity of a metal. The 
expectation value of the current in presence of a static field E leads always to zero 
conductivity because at thermal equilibrium no net current can flow. Indeed we have to 



take the appropriate limit, with a time dependent field and take CO ^ after. This leads 
to the generally accepted Kubo formula for the conductivity. 

In the superfluidity experiment, on the other hand, we have to consider the physical 
case when the relaxation time for the current becomes macroscopically large (infinite 
for infinite size) and the limit CO ^ after the limit L — > °° at finite temperature T leads 
to: 

" p 

Jdt<J(t)J(0) > 
o 



< J >r 



(l-p s /p)<K>v 



(12) 



where p s is the definition of superfluid density, being p the total density of bosons. 
Indeed whenever p s > only a fraction of the particles 1 — p s / p can be interpreted to 
have relaxed to the steady state in an experimentally accessible time. It is important 
to emphasize that this definition of superfluidity is experimentally testable but is not 
necessarily related to broken symmetry. Indeed in two dimensions, as well known, 
superfluidity can be detected[l 1] and on the other hand there is no long range order at 
any finite temperature. Whenever there is long range order the situation is indeed more 
conventional because the superfluid density is directly related to the helicity modulus of 
the order parameter []]]. 



FREE BOSONS 

In the free boson case, as shown in [1], it can be proved that p s coincides with the 
condensate fraction p c of particles that occupy the k = state with a macroscopic 
occupation. This calculation can be immediately generalized even for lattice models 
in any dimension D. Since the generalization is almost immediate we describe the basic 
steps in D = 1 for convenience of notations and write down the explicit expression of p s 
as a function of D only in the last equation. 

For free bosons the current J commute with the Hamiltonian and the linear response 
kernel a(0) is given by: 

a(0)=/3</ 2 >. (13) 

In fact the Hamiltonian, as well as J ©, are diagonal in k space H = Y,k £ k n k J = 
Y^k^o ^3k n ki where in the latter expression we have for convenience removed the k = 
vector in the summation because ^ = for k = (the derivative of an even function 
is odd in k). Following Ref.lll], in order to evaluate Eq. (fT3l) it is enough to compute the 
two body density matrix in momentum space: 

< n k n k , >= n B k 4 + 8k#[{4f + 4} d4) 

where 4 = l/[ e ^P( £ k — A*)) — 1] is me f ree boson occupation at finite temperature, and 
jl is the chemical potential used to require a given density p of bosons p = jY,k4 = 
p c + \ Y,k^o4> where p c is just the condensate density. In this way the evaluation of 
a(0) can be readily performed and simplified, by using that i) Y,k^o4l[t = ^ a S am 



because of the reflection symmetry ii) as noted in Ref-ffll [(nf ) 2 + nf] = — 1/ /3 

«(o)=-if^y^ as) 

k%\ dk ) d£ k 

We can now take the appropriate L — > °° limit to compute p s and p by replacing the 
summation l/L D Y,k^o ~^ I (ifp anc ^ obtain a closed form expression for p s (a simple 
integration by part is also left to the reader): 

pjp = (16) 

F / F f dkD » 

PcH 7T e 



d £ k 
dkl 



k x =0 



It is interesting that p s ^ p c in this case, but there is superfluid density only when there 
is a non zero condensate fraction and the other way around. 

Thus p s = for j5 < /3 C where f5 c is the inverse Bose-Einstein transition temperature 
that is finite in 3D but is infinite in ID and 2D. 



Zero temperature limit 

In principle we can take first the limit /3 — > °° for the kernel a at any fixed size L. 
As we have emphasized before, this limit cannot test superfluidity and indeed is related 
to another physical quantity, the zero temperature Drude weight as established by Kohn 
long time agoJH: ' " 

D c = lim ~ \im(<K> -a(0)) (17) 

L — >°° L, p — >oo 

that distinguishes a metal from an insulator, but not a metal from a superfluid. 

In the following the distinction between a Bose-metal from a Bose-superfluid is 
essentially analogous to the difference between a metal and a superconductor valid for 
electronic systems[13]: superconductors are obviously metal in the sense of infinite zero 
temperature conductivity but they also possess the non trivial property that the current 
can flow basically forever without dissipation at any finite temperature below T c . Clearly, 
within this definition, if T c = 0, we can speak about a Bose-metal in the ground state 
because there is no measurable superfluid density for any T > 0. 

In order to show that the limit /3 — > °° before the thermodynamic limit is incorrect 
for the detection of superfluidity, it is enough to realize that, for free bosons, in the limit 
/3 — > oo at fixed L the kernel Gf(0) — » because the current commute with the Hamiltonian 
and in the ground state J = so that cc(0) = /3 < J 2 > decays exponentially to zero for 
/3 — > oo due to the finite size gap 1/L 2 of the first excited state with non zero current. 
Thus we obtain that the Drude weight for free bosons is always finite and equal to < K >. 
Thus D c ^ in ID and 2D even though at any finite temperature 



<^>= a(0) 



for L — ► oo 5 implying p s = for any T > 0. 

In our definition therefore, a free ID or 2D Bose gas, is not a superfiuid but a Bose 
metal. This Bose metal however is too much idealized to be considered realistic because 
interaction is always present and it is known that an arbitrary small interaction changes 
the spectrum of the excitations from quadratic to linear in momentum, and condition iii) 
for superfluidity is not satisfied. Thus the issue of the present paper on whether Bose 
metals can exist in a stable phase is not solved by the free boson example. Bose metals 
should exist in nature only if a small physical perturbation of the Hamiltonian does not 
change the qualitative features of the unperturbed phase. 

Other definitions are known in the literature for a Bose-metal lfbUl . but appear much 
more restrictive definitions than the present one. 



MODELS WITH INTERACTION IN ID 

It is fortunate that the problem has been already studied in ID and we can use convincing 
results obtained in other contexts lfl5l H7I1 . The calculation of p s was done with the 
correct order of limits in id. In this work the authors considered ld-spinless fermions 
at half filling with nearest neighbor hopping, nearest V and next-nearest neighbor W 
repulsive interactions. This model, as well known, is equivalent to hard-core bosons 
with the same interaction coupling constants, because in ID hard core bosons with 
nearest neighbor hopping are simply related to spinless fermion models with the same 
current-current response functions. The spinless fermion model, in the gapless phase, 
is relevant for our discussion and it was clearly found that indeed p s > for T > 
as long as W = and V > 0. However the authors claim that, an infinitesimal small 
coupling W provides a vanishing p s at any finite temperature because-they argue- the 
model is no longer integrable by Bethe ansatz. In this case the zero temperature limit 
of p s does not coincide with the zero temperature Drude weight, that is generally finite 
in any gapless ID spinless fermion phase, because it is related to the low-energy zero- 
temperature properties of the model. 



If we agree with the conclusions of Ref.[ 15], that are based on calculations on periodic 



rings with ~ 20 sites, the model containing nearest and next-nearest neighbor repulsion 
is a Bose metal for any non zero V, W > 0. 

Indeed the conclusion of the work lfl5ll is more general and , translated in the boson 
language, implies quite generally that ID Bose metals do exist in the gapless phase. 
According to the authors conjecture, that is still under debate, p s = for T > in 
all models that are non integrable with Bethe ansatz (e.g. also the celebrated Bose- 
Hubbard model falls in this class if we extend this conjecture also to bosonic models). 
A more clear argument was given in Ref J16ll . where the absence of a finite Drude 
weight (p s ) at finite temperature was predicted in all ID models that do not have some 
conserved current. Essentially, in lattice models, the current can decay due to Umklapp- 
processes and a finite conductivity is expected at finite temperature, a condition that 
is incompatible with a finite p s (which implies a 8— function 00 = response and 
therefore an infinite finite temperature conductivity). The condition of integrability may 
instead allow for some conserved current, but it is also possible in principle that some 
conserved current can be realized even in non-integrable models. [18] Recently a more 



clear numerical evidence was also given that in a generic ID model with frustration p s 
is zero at finite temperature even in the gapless phase il7ll 

We do not want to enter in this subtle discussion on what is the right criterion that 
allows a finite p s at finite temperature, but from what is known so far, it appears that 
only very particular lattice models obtained with fine tuning of coupling constants can 
represent ID Bose superfluid and that the generic gapless phase is instead a Bose metal, 
at least for hard core boson models. Moreover, in this case, the superfluid phase obtained 
at particular coupling strengths do not certainly satisfy property (iii) and superfluidity 
may be detected only for suitable and very particular environment-system coupling. 



CONCLUSION 

We have formulated a consistent definition of superfluidity valid for lattice and contin- 
uous models in any dimensionality that relates p. s -the superfluid density-to the linear 
response current-current correlation calculated at finite temperature. This formulation 
agrees with the Pollock-Ceperley[l] one based on the winding number, provided the 
correct order of limit is taken: first the thermodynamic limit and then the zero temper- 
ature limit, relevant for ground state properties. In the opposite order of limits we have 
shown that the so called zero temperature Drude weight is obtained, but this can be fi- 
nite both for a Bose-metal and for a Bose-superfluid. The discrimination between the 
two can be obtained at finite temperature within the present formulation or by using the 



Scalapino-criterion that can be worked out directly at T = O fl 1 3fl . Both criteria coincide 
in the T — > limit for model systems where the solution is known, but the latter one 
cannot be applied in ID because it is not possible to define a transverse field in this case. 

The main conclusion of our approach is that in ID a generic gapless phase may be 
metallic and not superfluid, namely a very peculiar and interesting interacting phase - 
the Bose metal- with finite zero temperature Drude weight[19] but no superfluid density. 

In many recent papers the possibility to have this type of Bose metal has not been 
considered yet, especially in 1D[@, HI IZl 01, where it has been usually assumed that 
the gapless phase is superfluid. This attribute was originally used to characterize the 
classical 2D phase corresponding to the ID zero temperature quantum model. This was 
certainly correct but may be clearly misleading, because the superfluidity of the 2D 
classical model may be not related to the superfluidity of the corresponding quantum 
model at low temperature. 

In this work we have shown that ID hard-core boson interacting-sy stems should be 
Bose metals in the generic gapless phase, simply because for these models superfluidity 
cannot be detected at any T > (apart for the mentioned exceptions), even when the 
Drude weight is non zero in the ground state. 

Based on the above results, it appears possible that this Bose metal phase can be 
extended also to some model without the hard core constraint, because this constraint 
should not play an important role at low energy [7]. The seminal work by Fisher and 
coworkers on the mapping of the ID zero temperature Bose-Hubbard model to a clas- 
sical 2D model at finite temperature is perfectly valid as far as the critical behavior at 
the transition is concerned. However, in this mapping, the superfluid density of the clas- 
sical model (that can be finite below the Kosterlitz-Thouless transition temperature) is 



related to the Drude- weight of the quantum zero temperature model and not-obviously- 
to its finite temperature superfluid density. This quantity can be in principle different 
from the Drude weight, even at arbitrary small temperature, whenever the system is in- 
deed metallic and not superfluid. It is also clear that the analytical calculation of the 
"superfluid density" reported in Ref.([7]) for Luttinger liquids refers instead to the zero 
temperature Drude weight which is obviously finite, but does not necessarily imply su- 
perfluidity. On the other hand in the numerical calculation reported in Ref.[2], no finite 
size scaling is attempted at fixed temperature. Based on these considerations it appears 
important to improve further the numerical results of the ID bosonic models with soft or 
hard core constraint by using recent more accurate and powerful techniques [5], that can 
be extended to much larger system sizes. This may allow to establish more accurately 
the nature of the gapless phases of ID Bose models. 

In 2D close to a metal-insulator transition we have recently speculated[20] on the 
possibility to have a non Fermi liquid phase before the Mott-insulating phase. In the 
boson language this possibility can be realized whenever the phonon velocity c in the 
superfluid phase goes to zero before the Mott transition. In such a case an anomalous 
phase with finite zero temperature Drude weight but no superfluid density should appear 
between the Mott insulator and the superfluid. In this phase it can be also shown that 
there is no condensate, using a known relation based on the generalized indetermination 
principle. lElll In the language of spin liquids the Bose-metal is just a gapless spin-liquid 
of the type stabilized in the frustrated 7i — J2 model 12211 . Although in dimension higher 
than one all these examples are clearly not well established because they are based on the 
variational approximation, we believe that, since in ID the Bose (spin) liquid is stable at 
least in hard-core boson models, it is worth to consider this phase as a possible phase of 
matter even in higher dimensionality and especially in 2D. 
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